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ABSTRACT 


A model based on Kitaigorodsky's application of similarity theory 
and modified by McDonnell to forecast the mixed-layer depth was studied. 
The model applies during the warming season and is based on the theory 
of similarity. The parameters involved in the model were determined 
from bathythermograph data recorded at Ocean Weather Stations November 
(latitude 30N, longitude 140W) and Bravo (latitude 56 30N, longitude 
51W). Parameters were evaluated daily and grouped by months. Both 
seasonal and transitional MLD situations were treated. 

From these parameters, the form of the dimensionless function P(N), 
claimed by Kitaigorodsky to be universal, was determined by least squares 
fit to be best approximated by a second order polynomial. Forecasting 
equations involving P(N) were developed for each month and tested with 
data from the following years for both OWS ships. 

There is general agreement between the observed MLD and that found 
from the prediction equation based on the last year's P(N) for the same 
month and location. Month-to-month and spatial differences in P(N) cast 
considerable doubt on its universality, at least as determined by the 


parameters as currently defined. 
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}: Introduction. 

Extensive studies have been made on the ensonified bands of water 
in the sea in an effort to utilize better their potential for sound pro- 
pagation. Sound transmission in the upper layers of the ocean is for the 
most part determined by the vertical temperature regime. The need for 
more information about this thermal structure to increase the effective- 
ness of our ASW equipment and perhaps develop new ideas from this know- 
ledge is urgent. 

Various methods have been devised for forecasting the ocean thermal 
structure. Statistical predictions of the thermocline depth and sub- 
surface thermal structure have been the recent trend. The tools of this 
statistical approach have been either multiple linear-regression techni- 
ques or harmonic analysis of temperature cycles at various depths. 

The bulk of applied research, however, is still based on either 
dynamical models or on parametric empirical relationships. Inherent in 
dynamical analysis is the problem of mathematical complexity if all pro- 
cesses are considered; if simplifying assumptions are made, the reality 
of the model becomes questionable. Forecasting techniques based on 
empirical relationships are only locally valid with monthly or seasonal 
adjustments required. 

As pointed out by McDonnell [5] in his paper "Application of Simi- 
larity Theory to Forecasting the Mixed-Layer Depth of the Ocean", the 
theory of similarity represents an alternative approach in building a 
forecasting model. Kitaigorodsky [4] was the first to investigate the 
application of similarity theory as proposed by Monin and Obukhov [6] 
to predict the thermal structure in the upper layer of the ocean. In 


ie 


the development of this model, Kitaigorodsky assumed that purely thermal 
convection due to unstable density stratification was negligible and that 
vertical gradients of salinity are equal to zero. This imposed a season- 
al limitation on the resulting equations. Generally speaking, a stable 
density stratification exists in the upper layer during the warming season 
when the thickness of the nearly isothermal layer can be considered main- 
ly a function of wind mixing. Heat fluxes across the air-sea interface 
during the summer are positive (inward) and tend to build and strengthen 
the seasonal thermocline. 

With these assumptions, McDonnell applied the method of Kitaigorodsky, 
with some modification of parameters to develop a practicable forecasting 
model. In McDonnell's conclusion a recommendation was made that future 
research be applied in determining the form of the dimensionless func- 
tion P(N), inherent in the application of similarity theory, for various 
oceanic locations in order to test Kitaigorodsky's contention that P(N) 
is a universal function. 

The present author studied two distinct geographical areas using 
McDonnell's mixed-layer depth forecasting model in an effort to estab- 
lish the form of P(N). In this way, the form of the function P(N) could 
be better fixed and the possibility of its universality tested. Further- 
more, the practicability of McDonnell's model and parameters could be 
tested if realistic mixed-layer depths could be forecast using his 


method. 
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2. Review of McDonnell's model. 

McDonnell used data recorded at OWS Papa and the theory of similar- 
ity to develop a method of forecasting the mixed-layer depths associated 
with transitional and seasonal thermoclines during the warming season. 

The mixed-layer depth (MLD) was defined as the depth at which water 
first became 1C colder than the water at the surface. Usually, this 
depth could be accepted as the top of the seasonal thermocline. Transi- 
tional thermoclines were identified as those having a temperature differ- 
ence from the surface of less than 1C with a certain degree of permanence 
so as not to involve those of diurnal period. McDonnell considered the 
term ''MLD'' and depth of the thermocline synonymous and refers only to 
mixed-layer depths associated with either transitional or seasonal thermo- 
clines. Only secular, non-advective, and non-divergent processes were 
considered as influencing the MLD. Other processes contribute to MLD 
behavior which deviates from the model. 


The relationships developed by McDonnell are: 


MLD = Pin) ~~ = 


N= Q42 * 


where: = total heat present or excess heat in the upper 
wind mixed-layer, 
= representative maximum wind, 


coriolis parameter times 104 (2W@) sin@ x 104) 


= coefficient of thermal expansion, 


@es .S 


13 


= = a dimensionless function of N with the form 
of a first degree polynomial. 

To specify the form of P(N), equations (1) and (2) were solved for 
P(N) and N respectively. Then measured values of the parameters provid- 
ed 200 paired values of P(N) and N which were plotted together. The form 
of P(N) was found by curve fitting to this plot. Seasonal and transition- 
al MLD's were separately treated, a linear function P(N) being determined 
for each of these situations. 

McDonnell pointed out that, if the parameters chosen truly represent 
the controlling processes, then the plot of P versus N would have little 
scatter. Large scatter indicates assumptions were inadequate, e.g., 
divergence and advection are certainly important during some intervals. 

McDonnell's final equations incorporating the linear relationship 


for P(N) were: 


MLD = 29 = 25x10" wh 6) 
where 
P(N) = 2.9N = asxio"* as 
for transitional MLD and 
MLD = oa = 100 ar, (5) 
where 
Pin) = 3.89N - wx (07 (6) 


for the seasonal MLD. 


14 


i Area study selection. 

Several basic considerations governed the choice of the data used 
in this study. The first requirement was dependability, i.e., the 
measurements must be of acknowledged accuracy and recorded at a fixed 
location with appropriate frequency as nearly continuous as possible 
during the periods of interest; the second requirement was immediate 
availability, an important matter because of the limited time available 
for preparation of the study; the third requirement was that data be 
suitable to measure the phenomena the thesis attempts to describe, which 
means mainly that the effects of extraneous processes, such as internal 
wave activity, convection and advection be minimized or, at least, evalu- 
ated; and a fourth consideration was that the data come from geographi- 
cally and climatologically dissimilar areas and from different times so 
that the possibility of a universal function and its application to fore- 
casting could be examined. 

The requirements having to do with quality, frequency and continuity 
are satisfactorily met by the data from OWS ships; in fact there are 
few other sources for suitable data. The particular weather ships from 
which data were used were chosen in large part because of their being 
on hand in large quantities, thus providing economy of both time and 
money. 

Specifically, data available for the study represented two dis- 
tinct geographical locations, one in the Atlantic (OWS Bravo 56 30N, 51W) 
and one in the Pacific (OWS November 30N 140W). In addition comparison 
was available with Hepemmenis work at OWS Papa (50N, 145W). 

According to Tully [8], OWS November is contained in the eastern 
extremity of the large Subtropic Region in which the mid-ocean flows 
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are zonal and the waters respond to surface processes. Advection of , 
thermal regimes are minimal since no major current system is present. 
The location coincides with the mean position of the permanent Pacific 
anticyclone for the summer months, but effects of convergence in deep- 
ening the MLD can be estimated from Fofonoff's [1] mass transport cal- 
culations. 

OWS Bravo, however, located in the eastern sector of the Labrador 
Sea does not possess these ideal conditions. Random advective influ- 
ences may be present due to meandering of adjacent current ee 
Additionally, monthly mean patterns of atmospheric circulation show the 
presence of a deep low over this location; therefore horizontal diver- 


gence can be expected in the upper layers. To some extent, as at OWS 


November, this effect can be estimated. 


the West Greenland Current (warm) on the north and Labrador 
Current (cold) to the south could provide advective influences. ) 


16 


4. Calculation of parameters. 

The start of the warming season is evidenced by the onset of the 
seasonal thermocline; it remains in effect until after the autumn equinox 
when the seasonal thermocline settles to lower depths by convection and 
decays. Data to cover this period were selected from the months June 
through October. 

To determine the parameter MLD, observed values of MLD were plotted 
against time for each month, MLD's being read directly from the BT trace. 
Plots were made with the time interval three hours, the normal spacing 
of BT observations aboard ocean weather stations (OWS) ships. Both 
seasonal and transitional MLD's were plotted from the six to eight BT's 
available per day. A smooth curve representing the top of the thermo- 
cline or actual MLD was then sketched connecting the plotted points. In 
this manner an observation time with a missing BT report could be assign- 
ed an interpolated MLD. 

A mean MLD was computed from the four plotted MLD's during each 
twelve-hour interval starting with midnight Greenwhich. If more than 
one interpolated MLD was contained in the averaging process, the inter- 
val was not accepted. By assessing the MLD in this manner, the ambient 
variations due to internal waves hopefully were reduced. 

To determine Q, a BT trace was selected from each 12—hour interval 
studied that best represented the mean seasonal (and transitional, if 
it existed) MLD for that interval. The value of the parameter Q was 
determined from this trace representing the total heat in the uppermost 
layer. A step-by-step procedure for determining the value of Q is 


explained in appendix I with appropriate illustrations. The technique 
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used by the author represents a modification of McDonnell's method. 

The parameter W (representative maximum wind) defined by McDonnell 
is an average of the five highest winds reported in a 24-hour  yalsYiga 
that precedes the 12-hour interval of interest by up to 72 hours. 

The values of B ,» the coefficient of thermal expansion, are list- 
ed in table 24 as given by Sverdrup [7]. The value of the parameter 8 
is selected by entering table 24 with the surface temperature of the 
representative BT for the 12-hour interval being studied and the appro- 
priate salinity. 

Table 1 is a breakdown by OWS ship and month of the nearly 1500 BT's 
which provided the data for determining 628 paired values of P and N 
subsequently used in evaluating the form of the function P(N). Of the 
‘cotal paired values, 473 represent seasonal and 155 represent transi- 
tional thermoclines. 

The following equations were used to obtain the paired values of 


P and N from the parameters calculated for each 12-hour interval. 


P= (MLO) QAR (1a) 


Ww * 


Ne QG@Sr 


WwW (2) 


Tables 2 through 12 give the values of the parameters and the corres- 
ponding paired values of P and N for each observation time. The only 
irregulatiry in this process was September 1960 at OWS Bravo where the 
available data represented only the first 10 and last 11 days of the 


month. During the 10 day segment missing, the surface temperature became 


1 
(Normally eight wind reports are available in a 24-hour interval) 
18 
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less by 3.5C and the MLD increased by over 30 meters, indicating that 
other processes than those considered in the model may be involved. 
Therefore the data for September were split into two segments and treat- 
ed separately. 

With this change of season, the heat fluxes across the air-sea inter- 
face, although not computed, may well be negligible. During the follow- 
ing month, October, (as the cooler continental air masses became more 
prominent) instability mixing due to density increases created by 
evaporation may influence the depth of this isothermal layer. The in- 
fluence of evaporation, not considered in this model, would be indicated 


by the scatter in the paired values of P and N. 
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TABLE . 2 


PARAMETERS USED TO DETERMINE VALUES OF P AND N 
FOR JUNE 1957 AT OWS NOVEMBER 


12.9 31.3 22a 


DATE . W Qo Ge emp, MEDye Ts P. N 
(KNOTS) (Kg cal/om®) (METERS) . (°c) “x 10% § 
06015748 13.09 6.2" 9ey 27.9 seu 21.1. 1.48 ..94 
060157. 12.6 4.8 27.3 2.8 81.22. -77° 
060257. 11.0 5.8 27.7 99.21.97 » 4.07 | 
060257. 10.8 4.8 26 .6 21.7 1.57 8? 
060357. 10.2 4.9 24.6 > me 1.71 F298 
060457. 10.2 5.8 25.0 23.372. 1118 
60557 a 10.2 a2 Ga. 760 20.3 en S.8 9 23.3u8 2.83 1.50 
060557. 9.0 «747 688 23.8 030.9 922.8 M3.34 £1.74 
OG0657Me 6.0 mS. 9MAE.17°6 25.5 UmeS.7 9022.2 203.0lu.1.40 
QGU65TEN 10.6: 084. 3b. 05) 2428 wl0.0 22.8 ™ 2.05 11.21 
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061057mme SoG) 48.GmEE 1.71 25.7 (1288 9 22.2 4.16 11.95 
O61057me 16.86 9, sae? SSM 27.0 11,8 23.2. 5.04 .2.24 
Goll57me 6.@) &.0m 91.7808 25.0) mell.7 @ 23.5 9.5.29 §2.31 
OG1ISTME 7.07.6 2.0408 23.7 110.18 22.4 . 5.14°2.14 
G61257m 6.Qn 10.am 9.005) 26.7 8:5). 23.4 11.96 3.67 
06135708 07.0» L1s4m] 2.68% 2790 § 9.19 23.6 9.54. 3.38 
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061557. 11.4 1.68 1 c6ue 922.5 
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061857. 13.8 11:6 25.9\ 27 1 2. F171 
061957. 14.8 11.9 . 31.15 22.2 2.49 1.63 
061957 . 14.6 11.7 27.8 22.5 2.25 1.63 
062057 14.2 8.9 29,2 21 Te T.84 pl. 26 
062057 . 14.8 12.3 28.2 71.8 2.28 “1.63 
062157 14.8 11.6 27.5 22.1 2.14 1.59 
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W Q, Qt MID, MD TS 
(KNOTS) (Kg cal/cm*) (METERS) (°C) 


PARAMETERS USED TO DETERMINE VALUES OF P AND N 
FOR JULY 1957 AT OWS NOVEMBER 


P N 
“x 104 * 
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PARAMETERS USED TO DETERMINE VALUES OF P AND N 
FOR SEPTEMBER 1957 AT OWS NOVEMBER 
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TABLES (Cont'd) 

DATE W e. z., MLD. MLD TS 

(KNOTS) (Kg cal/em”) (METERS) (aC) 
092457 6.0 2.45 27.4 24 
092557 7.0 20.0 1.36 61.0 29.9 23 
092557 5.0 23.6 1287 «67.1 33.5 24 
092657. 10.6 27.2 £2.17 54.9 24.4 24) 
092757. 16.4 29.3 67.1 25 
O72e757 8 «4617:6 22.4 54.9 25 
092857. 17.6 26.2 61.0 25 
092857. 17.6 #£26.2 61.6 25 
092857. 17.6 #£28.0 67.1 24 
093057. 11.6 #£=22.6 45.7 24 
093057 11.6 22.6 45.7 24 
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TABLE 6 


PARAMETERS USED TO DETERMINE VALUES OF P AND N 
FOR OCTOBER 1957 AT OWS NOVEMBER 


DATE W Q Q MLD MLD TS P N 

S t, S t s ,S8 

(KNOTS) (Kg cal/em’) (METERS ) (°C) x 10 

100157 10.6 19.3 30.2 24.6 8.13 3. 
100157 10.0 24.3 36.0 24.7 FSoeyZ2 5. 
100257 10.6 22.2 35.6 24.4 TOS 4. 
100257 13.6 25.8 35.01 24.4 7.67 4G. 
100357 13.8 21.0 36.6 Be 5 S55 3. 
100357 15.9 227.7 44.2 24.7 9.93 4. 
100457 14.6 23.5 42.2 24.3 7250 3. 
100457 16.2 26.4 48.5 24.8 7.0/7 3. 
100557 16.0 22.8 39.8 24.3 S250 3% 
100557 lite? 23.2 50.1 24.4 6.95 3. 
100657 15.3 19.8 34.0 24.4 a5 2. 
100657 i4ee2°. 23.6 39.4 24.3 Tale +Oe 
100757 Wea 22.2 40.9 Dae Gop Ze 
100757 21.0 24.1 44.9 24.2 3.84 2. 
100857 22.0 23.0 42.4 24.1 sa 2 
100857. -_- 21.8 22.3 40.4 24.3 2398) 2. 
100957 22.0 28.4 46.3 24.4 “aze 2. 
100957 218 «621.6 37.1 24.5 2-65.52; 
101057 17.6 22.4 43.8 24.3 4.97 2. 
101057 14.8 20.9 41.7 24.3 6523" 3: 
101157 14.8 20.5 38.0 24.2 5.58 32 
101157 Hees 22.4 42.3 24.1 6.78 3. 
101257 14.2 26.1 44.7 24.2 9.07 3. 
101257 13.8 26.4 45.2 2a 2 9783 4. 
101357 10.0 19.6 46.9 24.3 4.42 4. 
101357 Gi? 222 43.5 24.2 17.89 5. 
101457 11.2 20.4 40.6 24. 1..-10s55: 3. 
101457 12.0 21.4 41.8 23.9 o.45 - 3. 
101557 12.0 21.8 44.3 23.5 WOn20 3. 
101557 12.0 19.0 39.6 22.8 foal 3. 
101657 Hees 21.7 45.3 22.8 Wea G&: 
101657 10.8 21.3 ea 21.8 11.30 3. 
101757 9.2 20.8 Cpe 22.8 15.36 4. 
101757 8.6 20.2 39.8 22.4 16.05 4. 
101857 8.6 20.0 44.0 21.7. 17207 «6G: 
101857 9.4 21.6 44.0 22.2 15.89 4. 
101957 7.0 20.1 41.9 22.2 25538 5. 
101957 10,6 2222 43.1 22.2 2.58 4. 
102057 1422. 2223 39.8 22.9 6.51 3. 
102057 15.2 25:5 42.7 22.3 6.95 3. 
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DATE W £2 

(KNOTS) (Kg 
102157. 16.6 27 
102157. 18.0 24 
102257 19.8 23 
102257 19.8 23 
102357 19.8 25 
102357 19.4 26 
102457. 16.8 25 
102457. 11.6 24 
102557. 12.6 24 
102557. 12.6 + 23 
102657 12.6 21 
102657. 12.6 22 
102757. = 10.0 25 
102757 10.6 20 
102857. 10.6 25 
102857. 10.8 25 
102957 9.4 24 
102957 3.0 722 
103057 7.4 21 
103057. 12.0 23 
103157. 16.2 27 
Posts? 87.2 "27 


NVUNMNUODKUUNNEUUHEWAYWNOOWON 


os 


Salven) 


TABLE 6 (Cont'd) 
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DATE 


060160 
060160 
060260 
060260 
060460 
060460 
060560 
060560 


060760 ' 


060760 
060860 
. 060860 
060960 
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061060 
061060 
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23.2 
Zo, 2 
16.4 
M2 
tye 
ie 
17.2 
TVor2 
14.0 
i522 
15.4 
17.4 
17.4 
17.4 
17.4 
Zi 
2136 
Z1pn6 
21.6 
2026 
16.6 
Mieke; 72 
19.8 
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PARAMETERS USED TO DETERMINE VALUES OF P AND N 


JUNE 1960 AT OWS BRAVO 


Qs Qt. MLD, MLD 
(Ke cal/em*) (METERS) 
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TABLE 7 (Cont'd) 

MLD. MLD, TS 
(METERS) @C) 
12 oe 

29.0 16.8 5.6 
35.1 14.6 5.6 
33.5 Tie © Se 
35 «1 18.3 6.1 
31.4 6.1 
32.0 6.1 
36.0 5.6 
33.5 6.1 
25.9 6.1 
27.4 6.1 
16.8 6.2 
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DATE 


070160 
070160 
070260 
070260 
070360 
070360 
070460 
070460 
070560 
070560 
070660 
070660 
070760 
070760 
070860 
070860 
070960 
070960 
071060 
071060 
071160 
071160 
071260 
' 071260 
071360 
071360 
071460 
071460 
071560 
071560 
071660 
071660 
071760 
071760 
071860 
071860 
071960 
071960 
072060 
072060 
072160 
072160 
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14.0 
14.0 
14.0 
12.6 
122.2 
nL) 64 
ise 2 
i. 2 
12.4 
12.4 
14.4 
18.2 
BZ 
SZ 
17.0 
12.6 
13 26 
13.8 
13.8 
13.6 
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FOR JULY 1960 AT OWS BRAVO 


Q Qs 
(Kg cal/om*) 


11.26 
5). 29 
i309 
H3.92 
13.64 
20.13 
£3.75 


15.05 . 


14.88 
14.77 
14.72 
13.24 
1322 
15.72 
2 oO 
10.35 
13.96 
15.74 
12.47 
11.01 
15.56 
1S faa 
16.58 


16.44 
15.75 
14.60 
O29 
15.52 
16.43 
14.13 
19.45 
16.45 


13.14 
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14.66 
17.57 
15.14 
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MLD, MLD, 
(METERS ) 
55.8 27.3 
57.4 26.2 
61.2 24.6 
63.4 34.4 
67.8 38.3 
57.4 31.7 
68.4 31.7 
65.6 32.8 
60.1 16.4 
7b. 27.3 
73.8 35.5 
79.3 30.1 
71.1 24.6 
65.6 27.3 
60.1 30.1 
54.7 27.3 
73.8 26.2 
76.6 30.6 
65.6 23.0 
61.2 12.0 
79.3 16.4 
65.6 27.3 
S20mam 33-9 

20.8 

76.6 23.5 
67.3: 26.2 
75.5 21.9 
67.8 25.7 
Folegl 27.9 
67.3 2153 
82.0 24.6 
87.5 32.8 
65.6 37.2 
60.1 35.5 
59.1 32.8 
71.1 36.1 
76.6 44.8 
62.9 S787 
36.1 

68.4 45.9 
19.7 

68.4 44.8 
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PARAMETERS USED TO DEVERMINE VALUES OF P AND H 


x 104 , 
5.90 “le 22 
Ja «=6R G3 
Jvae 14642 

WO323 1.67 
Wie 3 269 
6-14 1216 

10-83. 1,72 
Vi.3> 1.86 
11.65 1.98 
12.57 1280 
9.64 1.55 

6235 1.20 
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6.49 1.48 

4.603. LAG 

7ol@.° 1435 

10.84 1.67 
12.68 1.88 
8.60 1.49 

7.09 T1.3t 

13.50 1.93 
9.42 1.70 

12.94, 1092 
12.49 1.99 
10.65 1.88 
13.34 1.95 
12,03 1:96 
13:35, 2.07 
14.22 2.26 
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TABLE 8 (Cont'd) 
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MLD. MLD 
(METERS) 
Fe | 
54.7 
71.1 
73.8 38 
61.2 35 
74.1 42 
71 38 
13 
73.8 49 
27 
33 
68.4 32 
iS egey 30 
19 
76.6 a1 
21 
76.6 20 
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DATE 


080160 
080160 
080260 
080260 
080460 
080560 
080560 
080660 
080660 
080760 
080760 
080860 
080860 
080960 
081060 
081060 
081160 
081160 
081260 
081260 
081360 
081360 
081460 
081460 
081560 
081560 
081660 
081660 
081760 
081760 
081860 
081960 
081960 
08 2060 
08 2060 
082160 
082160 
082260 
082360 
082460 
082460 
082560 
082560 
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TABLE 9 


PARAMETERS USED TO DETERMINE VALUES OF P AND N 
FOR AUGUST 1960 AT OWS BRAVO 


Q Q+, MLD MLD, 1S P N 
(KNOTS) (Kg cal/em®)  (MmrERS) “ (°c) “x 104 ® 


3.00 1.43 


6.80 a Sie 8.9 

9.06 18.3 9.3 26/Smn 1547 
4.67 1 Sa 8.9 “o) 62 
8.40 Qiks 9.2% 2e260) 1 a9 
9.46 22.9 9.20 17756) ) 07 
15.85 3045 9.3 3.91 1.79 
10.25 24.4 8.9 2-507 27 
6.22 18.3 9.2m) 1285 715008 
351 B0.5 10.1 680mm 2523 
14.42 1.45 29.0 9.18 .9.45 9.77 2.77 
nr. 7.2 25.9 9.5m Sesee 2,1 
22.96 3.96 44.2 12.2° 9.7 10.56 2.99 
10.02 24.4 om lls 257 cee 1.43 
14.92 29.0 10,4ymea 65 2.00 
tO. 50 25-9 11, es 10@ 1.49: 
3 02 24.4 11. ee? Oommen 3c 
8.66 18.3 id, 102 93 
2.35 2). or 1 Ga . 70mm) 
8.14 oni 1. Zee. 7Omee 1. 29 
10.82 18.3 11, So. Came 300 
SoS 15.8 11,9 2.4 eee 50 
2,30 19.8 te). 7 Wee. See 2 sa 
13.06) pieces 22-90) Ome) ,7997 3,32) 1785 
11.66 21.3 11.67 2.75 1.65 
15220 29.0 AW; | es Tee 
14.64 a7 ls i Gue, 39%ee 0S 
16.50 30.5\ 1. Jee 9. COmmeon sO 
20.376. 6.55)38.4 @elcmomee 11.7 911,05 3.27 
9.62 ales) Wey ae) La bys 
13.06 25.0 i 7a « 22022 
B2.50 24.4 il: 6amo.07 1.95 
20.58 29.0 i. Me2.36 «1,473 
19.05 27.4 tO. Jee. 02M 5G 
10.52 e225) 23-5 One!) Geen, O4 91 
12,90 (eee 32.0 es ie Le 2.77 9 1.40 
PLS 39.6 1 hee Somes 44 
12.00 18.3 ] ee. 349 2). 50 
17.30 26.8 iM. 1S 257 2.53 
15.28 ley it. 7.00.9 2.28 
1525 25.0 i ee GA 22 
1S, Sie 251.6 ie Hees O30) 2,26 
9.35 18.3 17 1.99 1.36 
12 G2 meme O20. 7 3 See o4 | 2,31 
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TABLE 9 (Cont'd) 


DATE W a oh MLD. MLD tS 
(KNOTS) (Kg cal/em“) (METERS) (°C) 
082660 11.8 19.60 36.6 11.7 
082660 11.8 22.62 .90 39.6 9.1 11.4 
082760 11.8 21.70 37.8 11.7 
082760 9.8 16.88 .67 30.5 10.7. 11.7 
082860 9.8 19.62 33.5 11.9 
082860 9.8 18.78 32.0 12.1 
082960 9.4 19.07 1.00 32.0 11.3 12.2 
082960 8.8 13.82 25.0 11.2 
083060 10.6 15.00 1.22 27.4 15.8 10.6 
083060 16.0 16.90 22.9 10.6 
9ea160 YO. 17.15 27.4 ie 
083160 19.4 16.00 26.5 10.1 
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DATE 


090160 
090160 
090260 
090360 
090360 
090460 
090460 
090560 
090560 
090660 
090660 
090760 
090760 
090860 
090860 
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092660 
092760 
092760 
092860 
092860 
092960 
092960 
093060 
093060 


(KNOTS) (Kg cal/emé) (METERS) (°c) 
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16.6 
19.2 
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18.0 
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15.0 
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2520 
25.0 
Asp A8 
28.4 
28.4 
30.0 
30.0 
22.8 
20.0 
20.0 
bse, 
21. 2 
22.6 


TABLE 10 


FOR SEPT. 01-09, 1960 AT OWS BRAVO 
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TABLE 11 


PARAMETERS USED TO DETERMINE VALUES OF P AND NN 
FOR SEPT. 19-30, 1960 
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PARAMETERS USED TO DETERMINE VALULS OF P AND N 
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TABLE 12 


PARAMETERS USED TO DETERMINE VALUES OF P AND N 
FOR OCTOBER 1960 AT OWS BRAVO 


DATE W Qs Ot MD, MD, %% P, ,Q 
(KNOTS) (Kg cal/cm2) (METERS) (Cc) x 10° 

100260 25.5 19.75 Cre be 1S eel? 
100360 27.0 19.05 58.8 ARON S3m 1.07 
100460 25.0 18.62 oo 5 emer o.76.°~«1.02 
100460 25,0 18.80 56 5 fe Gl ee 14 
100560 20.0 20.40 54.3 sao 4.63" PaO 
100560 16.8° 17.78 © 59.4 PeGEES Meco 
100660 18.2 15.25 51.2 Se 4.00 1 
100660 19.6 17.55 59.4 6.0 5.08 1.38 
100760 19.6 19.95 62.5 6. 1M O7 eI S4 
100760 19.6 15.32 47.2 SASS, 14 O07 
100860 25.0 18.45 59.4 5:0 2.93 1,01 
100860 25.0 17.72 57.9 ey, 
100960 23.0 11.20 42.7 Se a dS meee 
100960 20.0 15.71 50.9 6.3.68 92.19 
101060 17.0 15.31 ow 6AOMNSS «411.36 
101160 20.0 10.38 39.6 od)” Sey 
101260 20.0 10.52 Bi 5 See 1.81 0.72 
101360 20.0 15.80 48.8 pele 2. 5SSuelead 
101460 22.0 15.55 47.9 5 Oe 2.57 OF 
101560 25.0 15.95 60.0 a eT ey 
101660 25.0 17.95 58.8 5.8 /2.82 99 
101660 25.0 16.50 Eaue nw ee se 
101760 17.0 18.10 62.5 home 6.520 10ae 
101760 17.0 . 13.26 51.8, poate Te MAC 
101860 17.0 16.62 62.5) 4 ll 6 ll 
101960 12.0 17.84 _ 62.5 \ 5 6M 12.93. On 
102060 20.2 14.71 aie ome 3 60m 1. 10 
102160 20.2 17.35 64.0 baleen S<00 1.30 
102260 20.2 17.24 58.5 ee 3 1.17 
102360 18.0 17.65 59.4 — ee 
102460 1.80 16.60 Same a Sk eee 
102460 Wy u7-0"weis..98 a Cuan Sesaene 138 
102560 13.0 17.30 61.0 5.7mm0.43.«1.83 
102660 20.0 18.52 67.1 ; Te MG 
102860 25.0 15.62 59.7 6 75 99 
102960 22.0 14.40: a) 6 gl Sableeme 98 
103060 22.0 14.57 62.8 amme | egies 
10326027600 w16.10 a? oi MeCeGGe Be 0 
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oe The form of the function P(N). 

A least squares computer program was used to determine the poly- 
nomial of degree K which best fits (in the least squares sense) M data 
points. The best fit among those polynomials tested (through third order) 
was for K = 2 for each of three groups of points representing about one- 
fourth of all paired values of P and N. The coefficients of the poly- 
nomial were then computed for each month and tabulated in table 13, 


P(N) having the form below, 

) ao 2 ¥ d 

P(N) = aN 4 AN 42, (8) 
The corresponding forecasting equation is 


we 
MLD = a,4Q soy * Dora gt (9) 


McDonnell's criteria for acceptable data limited the number of his 
paired values to only 22 pairs for transitional MLD's and 29 pairs for 
seasonal MLD's. These data, as a result, were from various months of 
the warming season during the years 1958 through 1962. Because of the 
small number of paired values and the grouping of the seasonal and transi- 
tional paired values, only a linear regression separately done for the 
two categories was justified. These are equations (4) and (6) of 
McDonnell; they do not necessarily represent the most likely form of 
the function P(N). 

The present author used both seasonal and transitional paired values 
together to obtain a single form for P(N). This was then incorporated 
into McDonnell's basic equation (1) and used to forecast both seasonal 


and transitional MLD's. Graph No. 1 represents the form of P(N) using 
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TABLE 


13 


COEFFICIENTS FOR EACH MONTH USED IN THE FORECASTING EQUATION 


(1-10) 


» (19-30) 


lise 


OWS November (June through September) 


OWS Bravo (June through September) 


»543 


-996 


OWS November 


OWS Bravo 
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1. 


-4.83 


1.289 


1.815 


-.228 


.023 


a second-order polynomial as the best fit for the paired values deter- 
mined by McDonnell at OWS Papa. 

Graphs No. 2 through 12 are the curves of the function P(N) as 
determined for each month. All paired values are plotted on each scatter 
dilegmem. 

The scatter of the paired values is relatively small for most months 
indicating that McDonnell's model may well contain the correct combina- 
tion of parameters. Usually the paired values of P and N for transition- 
al situations were found near the origin with little Pon During 
low wind conditions, the computation of P is very sensitive to small 
errors in wind speed which accounts for much of the excess scatter at 
large P. Additional scatter probably results from random fluctuations 
not removed by the averaging procedures SPienihed in section 4. 

One can see that the monthly best fit curves have a variety of 
slopes apparently indicating the non-universality of P(N). However, 
systematic deviations due to contaminating influences (e.g. divergence), 
but included in the computation of the paired values, may account for 
the variations in slope of each monthly function. By analyzing incre- 
mental changes in P and N associated with small increases in Q and MLD, 
general conclusions concerning the influence of divergence and advection 
on the paired values can be made. This analysis indicates that reduc- 
tion of the MLD by divergence or advection tends to diminish the slope 


dP/dN and vice versa. 


a (ezach No. 5 for September 1957 had 10 points which fell outside 
the scale. Graph No. 4 for August had one such point.) 
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Divergence of the Ekman transport was computed from the monthly 
Ekman transport at grid points in the vicinity of each location during 
the year studied. Meridional and zonal components of Ekman transport 
calculated by Fofonoff and Ross [1,2], were used for this. At OWS Bravo, 
maximum divergence was during August which has the least slope of any 
function for that OWS ship. The same correlations were noted at OWS 
November except that the divergence was negative.” 

Systematic deviations in the paired values as a result of advection 


could not be evaluated as easily. 


* (aad and August at OWS November were anamolous months in this 
respect.) 
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oy A possible universal function. 

The concept of a universal function P(N) as proposed originally 
by Kitaigorodsky was investigated by combining all of the 504 paired 
values of P and N for the months of June through September for both OWS 
ships. By the least squares best fit method, the second order poly- 


nomial for P(N) was found to be 
-4 
P(N) = 422K (ONt 4 2.25N — eB xI0 (10) 


with the resulting universal forecasting equation, 


Bes ety 
MLD = 422x108 Q41. 26 — Iba x10 sant an 


Graph No. 13 represents the function P(N), equation (10), with 
upper and lower bounds of one standard deviation of the residues. The 
residues are defined as the difference between the computed and original 
ordinates and can be interpreted statistically as the standard error of 
estimate of P. 

Graph No. 13 Hisc indicates the least-squares best-fit polynomial 
for each OWS ship during the same months June through September. The 
function P(N) from OWS November remains inside the statistical bounds 
indicating that the proposed universal function may be appropriate for 
that location. OWS Bravo, located in a more dynamic area, has a func- 
tion which exceeds the statistical bounds for high values of N. Pro- 
cesses not included in the model may explain this deviation. 

The function P(N) for each OWS ship is estimated from the data of 
only one warming season and may well be unrepresentative. Investiga- 


tion of other years may reveal a closer correlation between different 


By 


locations and times which would strengthen the idea of a universal 


function as well as improve the estimates of the constants involved. 
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7. Procedure for forecasting and testing. 


Equation (9) can be used to forecast MLD's over any length of time 
for which the parameters can be accurately predicted. Data such as 
were used to determine the coefficients in (9) were available for the 
following years at both OWS ships. A continuous day-to-day forecast was 
used to test the appropriate monthly coefficients for equation (9).7 In 
essence the forecast was a test of whether the curves P(N) for a given 
year and month were useful in predicting MLD's for the same month in 
some other year. 

All BI‘s available for the preceding 24-hour period were used to 
calculate a mean observed MLD. 7 The parameters 3 » Q, and W were con- 
puted by the same methods used in determining the paired values P and N. 
Using the parameters » Q, and W in the forecasting equation (9), with 
the proper coefficients for the month and location under study, a daily 
MLD was computed and compared to the 24-hour mean observed MLD. This 
process was continued day by day from the available data with the results 
listed in tables 14 through 22. A total of 169 forecasts were made, 20 
representing MLD. and 149 representing MLD .- 

Although forecasts for periods greater than 24 hours were not at- 
tempted, equation (9) is assumed to possess this utility. In an extend- 


ed forecast, a mean value representing the heat flux across the air-sea 


ee a small number of observations was available for June and 
July 1958 at OWS November. August data for the same location were 
missing.) 


nee comparison with the computed daily MLD, a 24-hour interval 
was necessary to provide additional BT data for averaging out non- 
periodic influences. ) 


a5 


interface per day could be applied to modify the parameter Q for heat 
accretion during the forecast interval. Monthly climatological data 
(Kimball [3] )are available for certain oceanic areas that list the aver- 
age net heat Elux per day. More important, however, is an accurate wind 
prediction. Its importance can be seen by analyzing the terms with the 
coefficients ay and ay of equation (9) from table 13, and noting the ex- 
pected changes in the parameters Q and W respectively. The average 
change in Q as a result of heat flux is at most about ten percent in a 
single day, based on approximately .4 Kg. enienc per day influx at OWS 
November, while the change in W may range from 0 to 30 knots during the 
same interval. When considering forecast changes in the seasonal MLD, 


the term involving the coefficient a, then becomes negligible. 


2 

Therefore, daily increases in Q were not considered essential in 
forecasting seasonal MLD's. The fact that wind through mechanical mix- 
ing during the warming season is usually the dominant factor in fore- 
casting changes of the seasonal MLD is clearly seen - assuming fluctua- 
tions created by internal waves have been averaged out. 


The possible universal function derived from all paired values for 


June through September was not tested by forecasting. 
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Deane W 
(KNOTS) 
062658 12.2 
062758 16.8 
062858 19.6 
062958 20.8 
063058 2O06Z 
Forecast 
Forecast 
071058 18.6 
071158 18,6 
OA a Ae 
O71355 12.8 
O7ss 10.6 
Owiess 10,0 


FORECAST OF MLD's FOR JUNE 1958 AT OWS NOVELNSER 


Qs Qt TS 
(Ke cal/em*) (°C) 


Liao 
1h. oe 
AAT (55 
LSP Mae 
ese 


20.0 
20.0 
eo .0 
20.0 
20.0 


TABLE 14 


VORECAST 
NLD, MLDy 
(J LTERS ) 
40.5 
39.6 
Leek 
Adel 
A762 


OBSE 


MLD. 


(MEGERS) ” 


48.8 
39.6 


43-5. 


4509 
46.6 


seasonal MLD's within one standard deviation 


rorecast+ 


DIFF 
oo 


DIFF 


(LTERS) 


~3.5 
30 
gla 
=. 

fe) 


(3.1 meters) 


seasonal MLD's within two standard deviations (6.2 meters) 


FORECAST OF MLD's FOR oULY 1958 AT OWS NOVEMPER 


9.26 
apreya 
9.20 
8.20 
9.60 
OmO2 


20.0 


TABLE 15 


38.9 
be.) 
3504 
32eL 
3501 
3569 


vi 
7 


32.0 
B05 
3807 
Sede 
37.2 


Le 


Forecast seasonal MLD's within one standard deviation (3.7 meters) 67% 


Forecast seasonal MLD's within two standard deviations (7.2 meters) 100% 


1 


(Negative values indicate forecast MLD's were too shallow) 
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DATE 


0901558 
090258 
090358 
090458 
090558 
090658 
090758 
090858 
090958 
091058 
091258 
091356 
O91458 
091556 
091658 
C91L758 
091858 
091958 
092058 
092158 
092258 
092358 
092458 
092558 
092658 
092758 
092858 
092958 
093058 


Forecast seasonal MLD's within one standard deviation ( 5.8 meters) 


TABLIS 16 


FORECAST OF MLD's FOR SEPTEMBER 1958 AT OWS NOVEMPER 


W 
(KNOTS) 


f- 
e 


} 
MUNN ®90OO 


° e e be ® e 


e 86@ & 
ODN OPANAODANNADAADRODIDOIUDOIOIOINDIOA 


SIRNA ROM DID 
oeeee 


e 


Beep ete Bee meet ek 


@ 


Q 


Ae 


12.60 — 


duly, 68 
Me 
As '75 
16.16 
15.60 
16.10 
lie 
Lefps.0 
tee 
16.70 
1690 
aS 70 
15.80 
L8e42 
19.77 
us 72: 
oop) 
20. ‘70 
1/09 
R25 
22.88 
2ST 
22650 
PATS 
Freeh 
Lape 
24.00 


S Qi5 
(Kg cal/om*) (°C) 


065 
.68 


aS 


2303 
owe 
220 
2866 
2303 
Ee 
2303 
2309 
e369 
2303 
AONE: 
A030 
pee 
coil 
aes, 
Ame: 
ele 
ZOeL 
2 al 
Reel 
e229 
ae ol. 
Bee, 
2298 
Pere | 
re 
PIE 
Co 
ae | 


FORECAST ORSLRVED PORECAST 
MLDg LD, MLD, MLD, DIFF, DIFF 
(METERS) (METERS) (ETERS) 

£917 37.5 ere 
50.6 39.0 11.6 
ee © 36.6 $.0 
BO. 3 38.4 1.9 
Lae) ISS ae 
41.8 10.5 1.3 
32.6 LOe2 = ILS 
36.3 33.1 ~ 1.8 
34.8 10.0 36.6 Ge, 26158 9 
3085 Oe ST. ON) =" Llc so 
Dees Boao a 1.6 
Phos Wl 39.6 "7 
L3: SORE Rae 
LbeY 39.6 ae 
Loe We See Loa 
Api ‘7 Ae. 7 Lg 
51.8 icra 7.6 
50.6 £19 eal] 
ae. 46.6 Die 
5.6 ree Jie 
L5e5 ER 1.3 
50.1 vec? 1.9 
ow TE way ui haa 7 
57.9 Sees Cel, 
58.6 51.8 6.8 
5696 Sees ore 
SARS 5.9 oF 
52.5 51.5 a) 
Sie 5 py Oa a Bele 
72% 


Forecast seasonal MLD's within two standard deviations (11.6 meters) 97% 


58 


TE 17 


FORECAST OF MLD's FOR OCTOBER 1958 AT OWS NOVANUER 


FORECAST OPSERVED FORECAST 
DATE W QO. Qi. aS MD, MID, MD, MID, DIFF, DIFF, 
(KNOTS) (Kg cal/em*) (°c) (METERS) (ETERS ) (1V2cTERS) 

HOG-50) 1).6- 21.85 22.0, paleo 47.9 ~ 6.9 
100258 Sey eeepc = 22.9 =e Lee — 5.9 
TOGs50 7.6 23.90 Pa me TS Pt! 549 — 2e2 
100658 6.0m 20265 Zeet Mice 5 se — Let 
100755 6.05 22515 2365) lee SEAS. Za 
100858 7 Opis 2352. 6550 BOs Bed. 
100958 10.4 22.40 23.30 oe Siiee ~ §.3 
101058 13.8 26.35 P3608 Bee ey, ~ 5.0 
101158 $11.0 25.65 2301. ieee 5469 — 3.6 
101258 Spl | AGING ZIRE Ie Sinai 57.9 3.2 
101358 6.0 25.00 2c Doe ep. Aes) 
101458 6.0 28.52 Daeg (One Ss 49 
101558 1 2umeee 305 2266S ie 5208 = 2.2 
101658 Oe 25400 2260 Dene SNe) ~ 1.8 
101858 19.6 28.45 22.6 lee 60-4, -~ 9.1 
HONesoe) 15.2. 25610 geel dere a8) - 9.1 
OOS Loe 27 oe 22 Ome 57.9 - 6.7 
HOA Se Le 21g ke 2209 50.3 Slice - 7,0 
1022) yu) 2 O16 9D 2s Bae 61.0 —- 7.2 
1024538 26.0 25.93 22.2 7 toe wer 59. $13.4 
102555) 16.0.626,10 21.2 Lone Su SS) —11.9 
102655 15.2 26.10 2261 Valor 5409 - 8.6 
102758 15.2 28.60 ; 21.6.9 oles 64.0 ~12.2 
HOZS56.. 12.0 26562 22.0 50.5 62.5 ~12.0 
102958 16.8 28.4¢C 2 Soe 61.0 -~10.7 
103058 27.8 29.60 2eed Ween ey ial - 9,0 

AS Ae Gis PEE py one? 65.5 —l1.3 


LO315& 


Forecast seasonal MLD's within one standard deviation (3.2 meters) 17% 


Forecast. seasonal MLD's within two standerd deviations (6.4 neters) 43% 
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FORECAST OF MLD's FOR JUNE 1961 AT OWS BRAVO 


FORECAST OBSERVED FORECAST 
DATE W Q,  Q&, TS MD, MD, MID MID, DIFF, DIFF, 
(KNOTS) (Ke’cal/em*) (°c)  (kévers) (METERS) (METERS) 
06195, 26.0 6.35 6.1. wel 36.6 eo 
062061 23.0 3.79 0.7 Vee 3249 Ped 
062161 23.0 3.25 6.1 33.9 Bren 6.5 
062261 17.2 2.38 6.7 25.3 152 10.2 
062361 16.8 4.58 cya 2374 AS alte 
062464 16.4 4.08 6.7 28.0 BIAS acd 
52501) 9 1oscmeee.. 16 Oey 0.7 Oe =e 
062661 23.8 4.26 6.1 aoe 2 oye 5.0 
062761 23.8 3.56 come seis 26.2 9.4, 
062861 19.6 5.87 .48 6.1 33.4 14.8 33125 9.0 °°53.3 “Same 
062961 15.8 5.07 .59 6.7 29.1 15.8 “30:5 9.1 <4 Gee 
06306. 10.0 5.8 .68 6.7 23.0 12.3 29.6 91 -6.6 3.2 


Forecast seasonal MLD's within one standard deviation ( 6.6 meters) 53% 


Forecast seasonal MLD's within two standard deviations (13.2 meters) 100% 
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TABLE 19 


FORECAST OF MLD's FOR JULY 1961 AT OWS BRAVO 


FORECAST (ORSERVED) FORECAST 

DATE W Qs Qe TS MD, MLD, Mig HED, DIFF, DITF, 
(KNOTS) (Kg calfem*) (°C) — (KETERS) (METERS) (VirEns 

@70O1GL* 27.6 5.34 6.0 3527 3240 Leo 
070261 28.9 8.04 6. 3a 36.9 is © 
070461 28.0 12.38 6. Lae Fees 61.9 nl od, 
070561 20.0 10.54 6.198 5ee2 Coat URES 
@7@commm 15.2. Seon 6.2 49.3 AS A: 6.0 
070761 20.8 9.16 6.1 oe 48.8 5.2 
O70S61 20.8 11.32 6.1 Mela 50.0 ay) 
O7096L 18.2 9.46 6.7 ie o2ae L769 Lag 
O7106ie 12.0 - 16,48 7.2 66,2 5409 lés 
G71 Gime? sO. a2 copes 3205, 5.6  5A;Cm lo nemo) «1262 Pll es oer 
O71 2ClM 1 2,0melc,o8 2.7% 5.8 62.80 tercmmeoc.e 18.3 Pe ene ae 
@71361 5 18.0mells.25 3.10 5.6 °73.60) 15eemeoo™ 32.0 Syl eden yal 
G7 7,6 ee 2060 5.72 5.6 25 ee Bom") ZO 
O7156L 20.0 5.85 5.0 2662 Ary. Pee 
C7Cole 1520 Oreo 5.6 37.6 oleae 6.2 
C7 7Ole 1560 Sette Del Bheek BOS eh 
073061 8.0 S Gem 7.2 35.9 26.2 9.7 
073161 10.0 Coes 7.2 One Pee Be yyy 


Forecast seasonal MLD's within one standard deviation ( 9.9 meters) 75% 


Forecast seasonal MLD's within two standard deviations (19.8 meters) 100% 
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TABLE 20 


FORECAST OF MLD's FOR AUGUST 1961 AT OWS BRAVO. 


FORECAST OBSERVED ‘ FORECAST 
DATE W 2%  &, TS MD, MD, FID, MID, DIFF, DIFF, 
(KNOTS) (Kg cal/em“) (°C) (METERS) (METERS) °  (12:'TERS) 

0161 18.0 7.10 9.5 16.8 pre ~ 1.6 
CS036L 16.0 6.52 8. Sl? 22.9 - O92 
0S0Z61 176 6.83 ge Lee USI AY Lhek ~10.7 
080561 12.0 1.84 6.7 

080661 12.0 1.50 7.8 

O€076L 10.0 150mm 7.8 

080861 13.8 ie Com,S 

080961 17.4 3.14 35 MA pie 272.0 
OS1061L 20.0 Pan 7.6 

081262 20.0 8.14 67 eyo 25.0 GAS 
OSGi ieee. 1500 Gay ton 32.0 15.3 
CS1A61 14.4 7.95 Cy ee Plage ~13.3 
O81561 14.4 5.94 yee 2 PEGE mG 
OS1Geiy 11,0 . 6.775 pla dL Zloty —l0;2 
081761 12.0 8.31 7 Sa tee Pte = hs 
OS1861 19.2 17.65 1 cme Tee 20 - 7.3 
061961 19.6 8.15 S. Siege Oink e626 
062061 19.0 8.80 7.6 18.2 29.0 =1Ove 
082161 19.0 11.12 Bro ete 30.5 ElOag 
Oo22ole Secu On2e soe Ley 30.5 =n 
082361 15.8 12.14 Sy aera 32.0 - 9.9 
082461 19.8 16.20 SEN As: 36.0 - 6.7 
082561 23.6 15.32 8.3 28.1 Bez ~ 9.2 
082661 27.0 17.10 9.3 331 BS Al ~ 5.0 
OS2761 24.8 16.12 9.3 30.9 39.0 ems e 
082661 22.6 16.02 8.3 28.5 35.1 = 656 
082961 20.5 18.95 93D e2. CM La -~ 4.3 
083062 15.0 14.40 S26: Meiad: een 1028 
‘olen Ile aWanyie 8.9 eee 29.6 ~ 5.5 


Forecast seasonal MLD's within one standard deviation ( 5.6 meters) 17% 


Forecast seasonal MLD's within two standard deviations (11.2 meters) 83% 
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TABLE 21 


FORECAST OD MLD's FOR SEPTEMBER 1961 AT OWS BRAVO 


| | FORECAST OBSERVED FORECAST 

DATE W Qs Qt. TS MLD, MUD, MLD, MLD, DIFF, DIFF 

(KNOTS) (Kg cal/em) (°C) - (METERS) j (METERS) * (uirers) © 
090161 17.6 13.39 8.9 19.0 | 274, - 8.4 
090461 20.0 13.53 8.3 18.9 274 - 8.5 
090562 18.0 16.30 8.9 23.0 30.5 Te 
090661 17.2 18.00 9.4 26.6 32.9 - 6.3 
091061 19.0 18.20 708 23.8 38.1 -14.3 
091261 15.0 23.05 8.3 32.5 Pt, 2 e Sy, 
091361 15.0 25.80 853°936.7 (eee 25.2 - 8.4 
091761 15.0 22,60 8.4 31.5% ~ , Abed - 9.3 
091961 20.0 24.89 7e8 3269 ‘4567 -12.8 
092361 16.4 26.30 768 3545 $769 =23.4 
092461 19.4 26.60 8.3 37.1 5409 -17.8 
092561 20.0 ah,095 §.3 3he5 54-9 ~20.4 
092661 20.0 26.20 8.3 36.4 | 57,9 eine 
092761 15.0 26.30 $.3 37 oo 67el =-29.6 


Forecast seasonal MLD's within one standard deviation (11.5 meters) 40% 
Forecast seasonal MLD's within two standard deviations (23.0 meters) 80% 
TABLE 22 
FORECAST OF MLD's FOR OCTOBER 1961 AT OWS -BRAVO 


10026 28.0 24.85 me y.2 Ose 7Oue Zieee 
100361 28.0 23.65 72 ©9129 me 76.2 67 
100461 20.0 18187 6.7 659 Dons lets 
100561 18.0 16.42 568 5163 C27 8.6 
100661 19.0 17.21 6.7 60.2 51.2 9.0 
- 100761 22.0 17.40 5.6 55.0 53.9 ce 
100861 25.0 17.65 5.6 57.8 © 51.8 6.0 
101161 22.0 17.68 5.6 655.7 67.1 ~11.4 
101261 26.0 19.27 5.6 61.9 64.0 = 2.2 
101461 20.0 13.55 CASI yee 64.0 -18.7 
5.6 49.0 67.1 -18,1 


101561 20.0 .15.42 
Forecast seasonal MLD's within one standard deviation (10.6 meters) 45% 
Forecast seasonal MLD's within two standard deviations (21.2 meters) 100% 
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8. Evaluating the results. 

Table 23 is a condensation of the statistical analysis of pre- 
dicted MLD . in relation to the observed MLD |. Deviations of the fore- 
cast from the observed MLD are compared with the standard deviation (Q ) 
of the daily mean of the observed MLD. for each month. Statistics were 
not obtained for transient MLD situations since too few of these occurred 
during any month for a statistical analysis. Persistence forecasts from 
day to day were used for comparison. 

Except for the month of October’, OWS November had a large per- 
centage of forecasts (72%) within one © , which is significant in that 
the average CY (5 meters) is small. 

For the same months at OWS Bravo only 40 percent of the forecasts 
were within one @ (9 meters). The inability of equation (9) to fore- 
cast accurately the MLD may be related to factors, such as divergence, 
not included in the model. Use of additional paired values P and N for 
each month should improve forecasts based on the resulting function P(N). 
Extension of the monthly study into other years should bring about fur- 
ther improvement, as random contaminating processes are smoothed out by 


increase in sample size. 


1 October was omitted to avoid months containing possible convec- 
tive mixing.) 
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TABLE 23 . 
COMBINED STATISTICAL ANALYSIS OF FORECASTS FOR SEASONAL MLD's 


% FCSTS % FCSTS o 

' MONTH . YEAR . OWS # OF FOSTS WITHIN ONED WITIIN TWOO = (METERS) 
June 1958 + November 7 80 : (50) 80 (75) Bem 
July 1958 November 6 67. (60) 100 (80) 3.7 
Septe 1958 November 29 72 (100) 97 (100) 5.8 
Oct. 1958 November 30 17. = (62) 43 (83) Bee 

' June 1961 Bravo , 12. dee mane) 100 = (100) 6.6 
July 1961. Bravo 12 58 (82) 100 (82) 9.9 
Aug. 1961 Bravo 24, 17. (91) 83 (100) 5.6 
Sept. 1961 «Bravo © 20 40 (100) 80 (100) 12.5 
Oct. 1961 Bravo ta 45’ (85) 100 (90) 10.6 
Overall average of forecast seasonal MLD's within one © 45 (82) $ 
Overall average of forecast seasonal MLD's within twg © 81 (92) Z 


al , 
(Values in parentheses are statistical analysis of forecasts by persistence.) 
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9. Conclusions and acknowledgement. 

As a result of this study concerning the application of a pro- 
posed mixed-layer depth forecasting model, the following conclusions 
can be made. 

(1) Persistence gives the best short term prediction of MLD 
in the locations studied. If no recent observations are available, 
predictions utilizing a previous year's P(N) and accurate wind fore- 
casts are useful. 

(2) The dimensionless coefficient P(N), inherent in the ap- 
plication of similarity theory, is best approximated by a second-degree 
polynomial. 

(3) A single function can be used to represent P(N) for both 
seasonal and transitional MLD's. 

(4) During the warming season, changes in the MLD are mainly 
influenced by variations in the wind speed. 

(5) The concept of a universal function P(N) proposed by 
Kitaigorodsky may be valid, but its determination requires consider- 
able refining of existing data to remove contaminating influences. 

For his invaluable aid in the preparation of this manuscript, the 
author is deeply indebted to Associate Professor J. B. Wickham, 
Department of Meteorology and Oceanography, U. S. Naval Postgraduate 


School, Monterey. 
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APPENDIX I 


METHOD USED FOR DETERMINING THE PARAMETER Q 


The parameter Q is defined as 


Q = ety (AREA) » \O" 


-c 
where the factor (AREA) is given by the integral Eins zat, qT, and T, 
1 


being the temperatures of the "isothermal" layer (see fig. 1, slide 1) 
below and above the thermocline (either seasonal or transitional), and 
Z is the depth from the surface to the temperature curve. Density is 

represented by - and a is the specific heat at constant pressure. 

In evaluating the factor (AREA), the most difficult step is the 
choice of T,- It is that temperature, where the water becomes isother- 
mal or nearly so. The isothermal condition may continue to great depth 
or exist in only a thin layer between temperature gradients. Frequently 
this layer is difficult to distinguish, in which cake reference must be 
made to adjacent BI slides to establish at least a nearly isothermal 
condition. In any case the subjectivity in calculating Q by this pro- 
cedure probably contributes to scatter of the curves P(N). 

Once Ty, and T, are determined, (AREA) is found by replacing Jy, 2a 
by an equivalent rectangle with the area 7, (T,, - TT): The depth of Z% 
is determined by a horizontal line drawn through the thermocline such 
that equal areas will result above and below Z (see fig. 1, slide 4). 

For OWS November during the warming season pc, = .975 (enivcene) 


for an average Salinity of 32.5°/00 and can be considered constant. 


For OWS Bravo — = 1.01 (cal@un:) for an average salinity of 34.5°/oo. 
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A constant factor was calculated that included p°, and a change of 
dimensions (from British to Metric and from Fahrenheit to Centigrade) 
enabling direct computation of Q from the BT slide. This factor was 
1/6.05 for OWS November and 1/5.9 for OWS Bravo. 
A sample calculation of Q., from slide 4 follows: 
1. Determine the difference in temperatures between T, and 


2 
T,- (13.8°F) 
2. Read the depth of the horizontal line Z. (150 Ft) 
3. If this slide were from OWS Bravo data, divide the product 


of steps 1 and 2 by 5.9, giving 


Q, = (13.8)(150) x 107) Seesise Eteecaliena) 


Calculations of oa are done in the same manner and usually are an order 
of magnitude less than Q.: 

This method outlined represents a modification to McDonnell's 
technique. He constructed Ty so as to intersect the BT trace at 200 
meters (656 feet). This method soon became unreasonable in evaluating 
Q for two reasons. First, excess heat in the uppermost layer was poorly 
Popredtiited: Q represented the excess heat in the layer above 200 
meters. Secondly, Q could be evaluated realistically only on slides 
from deep BT's which are seldom used. The present author's method, al- 


though subjective, better represents the excess heat in the mixed-layer 


under study. 
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TABLE 24 
COEFFICIENT OF THERMAL EXPANSION (4x104) OF SEA WATER 
AT SEA LEVEL FOR DIFFERENT TEMPERATURES AND SALINITIES 


SALINITY 0/00 


TEMPERATURE (°C) 
jt 
& 


ti 


30 31 32 33 34 
iipgay 1.04 1.06 1.08 aig 
1.12 1.15 rs yeag 1.22 
2s 1.26 1.28 1.30 1.33 
1.34 1.37 1.39 se 1.44 
1.45 1.48 1.50 1.52 1.55 
1.57 1.59 1.61 1.63 1.65 
1.67 1.69 R72 1.73 1.75 
1,77 1.80 1.82 1.83 1.84 
1.87 1.89 1.91 1.93 1.94 
1.97 1.99 2.01 2.02 2.03 
2.06 2.08 2.09 Pysaiat 2613 
2.15 2.16 2.17 2.19, i 2021 
2423 2024 2.26 2.28 2.30 
ALY: 2.33 2.35 2.37 2.39 
2,41 a 2.42 2.44 2.46 2.47 
2.50 A 2.53 2.55 2.56 
2.58 2.59 2.61 2.63 2.64 
2.67 2.68 2.69 2.70 By 
2.75 2076 2618 2.79 2.80 
2.83 2.84 2.86 2.87 2.88 
2.92 2.93 2.95 2.96 
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